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Contrary to the recently published results by Daerden and Vanderzande (Phys. Rev. E 53, 4723 
(1996)), we show that the time correlation function in the random neighbor version of the Bak- 
Sneppen model can be well approximated by an exponential giving rise to a 1//^ power spectrum. 
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Recently, an exact solution of the random neighbor 
version of the Bak-Sneppen model was presented by de 
Boer and co-workers They derived a master equa- 
tion for the probability P„ {t) that n of out of N numbers 
have a value less than a fixed value A at (discrete) time 
t. In the limit N ^ oo and A = 1/2, P„ has the scal- 
ing form P„(i) = -^/(x = n/VN,T = t/N). Inserting 
this expression into the master equation gives the follow- 
ing Fokker-Planck equation for f(x,T) with a reflecting 
boundary at a; = 0: 
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Consequently, the random neighbor version of the Bak- 
Sneppen model for — > oo is just an Ornstein-Uhlenbeck 
process, i.e., Brownian motion in a parabolic potential. 
Given the initial condition /(x, 0) = 5{x~y), the solution 
is 
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FIG. 1. Plot of the correlation function given by Eq. |^. 

It follows for the autocorrelation function G(r) of the 
time signal x(t) 



G{t) 



1 



■ exp 



-a\r\ 



(3) 



where a = I. This directly gives the power spectral den- 
sity S{Cj) via a Fourier transform of G(t) 
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Going back to the unsealed variables leads to 
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which is only valid for low frequencies. Hence, the power 
spectral density of the signal n{t) decays as 1//^ and for 
very low frequencies it even becomes constant. However, 
the above calculation was carried out without applying 
the boundary condition at x = 0. Nevertheless, it is 
already clear from a physical point of view that the func- 
tional form of G(r) will not change drasticly by incor- 
porating a reflecting boundary. This is supported math- 
ematically by the fact that one simply has to use the 
method of images. This was done in Ref. ||^ giving 

G(r) = ^[1 - exp-2-]3/2[F(l, 2, 3/2, r_(r)) 

+F{1, 2, 3/2, r+(T)) - F{1,2, 5/2, r_(T))/3 
-P(l,2,5/2,r+(r))/3]-^, (6) 

where F{a, b, c, z) is the hypergeometric function and 
where r±(T) = ^[l ± exp(— r)]. In Fig. |^, G(t) is 
shown. We clearly flnd an exponential behavior with 
a = 0.869 ± 0.008 for 0.1 < r < 10 giving rise to a power 
spectral density as in Eq. (^ . This is conflrmed by a nu- 
merical Fourier transform of Eq. (^ (see Fig. |2|). Here, 
it has to be noted that G(t) is an even function, i.e., 
G(t) = G(— r). This also ensures that the Fourier trans- 
form 5 ((D) is real. 

In the case of the Bak-Sneppen model with one next 
neighbor, we also cannot conflrm the results presented 
in Ref. j^. A direct simulation of the time signal gives 
Sp{uj) oc 1/w^'^ over 2 decades for a system size of 
N = 8192. Hence, although the power spectral density 
in the Bak-Sneppen model decays as a power law, the 
exponent is far from one. This is also true for a different 
definition of the time signal 

In conclusion, there is no sign of 1// noise in the random 
neighbor version of the Bak-Sneppen model and even in 
the next neighbor version there is no l/f noise in the 
strict sense. 
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FIG. 2. Plot of the power spectrum given by a numerical 
Fourier transformation of Eq. ^ The solid line with exponent 
— 2 is drawn for reference. 
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